What is claimed is: 

1. A method for automatically segmenting lung 
5 nodules in a three-dimensional (3D) Computed Tomography 
(CT) volume dataset, comprising the steps of: 

receiving an input corresponding to a user-selected 
point near a boundary of a nodule; 

constructing a model of the nodule from the user- 
10 selected point, the model being a deformable circle having 
a set of parameters E that represent a shape of the nodule; 

estimating continuous parts of the boundary and 
discontinuities of the boundary until the set of parameters 
fi converges, using dynamic programming and Expectation 
15 Maximization (EM) ; and 

segmenting the nodule, based on estimates of the 
continuous parts of the boundary and the discontinuities of 
the boundary. 

20 2. The method of claim 1, wherein the set of 

parameters 6 = [O, s] T , O being a position of the model, s 
being a scale of the model, and T being a transpose of a 
vector corresponding to the position O and the scale s of 
the model . 
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3. The method of claim 2, further comprising the 
step of representing the boundary as a sum B, where 

B= (U,- B ci ) LKU, B dj ) , B ci represents continuous parts of the 
boundary and B d j represents discontinuities of the boundary. 

5 

4. The method of claim 3, wherein said estimating 
step includes the steps of; 

estimating the continuous parts B C i of the boundary 
based on a Maximum A-posteriori (MAP) estimate according to 
10 an equation B C i = arg maxB C i p{B C i%I, E) , I being a slice from 
the 3D CT volume dataset; 

estimating a MAP density as p{B c £l, E) = 1/z 
exp (E E (B C i) ) , E E (B C i) being a sum of internal shape and 
external image energies, and z being a normalization 
15 constant; 

minimizing the sum of internal shape and external 
image energies E & {B C i) using a time-delayed discrete dynamic 
programming method; 

connecting the continuous parts B c ± of the boundary to 
20 obtain an estimate of the discontinuities B d j of the 
boundary; and 
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updating the set of parameters £, based upon a circle 
fitting method being applied to the continuous parts B ci and 
the discontinuities B d j of the boundary. 

5 5. The method of claim 1, wherein said constructing 

step comprises the step of increasing a radius of the 
deformable circle until the radius contacts high gradient 
points in the 3D CT volume dataset . 

10 6. The method of claim 5, further comprising the 

step of pre-processing a region-of -interest that 
encompasses the user-selected point using an Expectation 
Maximization (EM) based method, so as to classify and 
remove a calcification from the region-of -interest . 

15 

7. The method of claim 6, wherein said pre- 
processing step removes the high gradient points that 
result from the calcification of the nodule. 

20 8. A method for automatically segmenting lung 

nodules in a three-dimensional (3D) Computed Tomography 
(CT) volume dataset, comprising the steps of: 

receiving an input corresponding to a user- selected 
point near a boundary of a nodule; 
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constructing a model of the nodule from the user- 
selected point, the model being a deforraable circle having 
a set of parameters £ that represent a shape of the nodule, 
where £ = [0, s] T , 0 being a position of the model, s being 
5 a scale of the model , and T being a transpose of a vector 
corresponding to the position 0 and the scale s of the 
model ; 

representing the boundary as a sum B, where 
B- (U,- B ci ) U(U,- B dJ ) , B C i represents continuous parts of the 

10 boundary and B d j represents discontinuities of the boundary; 

estimating the boundary, wherein said estimating step 
includes the steps of; 

estimating the continuous parts B C i of the 
boundary based on a Maximum A-posteriori (MAP) estimate 
15 according to an equation B C i = arg raaxB c i p(B c iZ,I, E) , I being 
a slice from the 3D CT volume dataset; 

estimating a MAP density as p(B C i^I, E) = 1/z 
exp (E s (B C i) ) t E E (B ci ) being a sum of internal shape and 
external image energies; 
20 minimizing the sum of internal shape and external 

image energies Eg(B C i) using a time-delayed discrete dynamic 
programming method; 
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connecting the continuous parts B C i of the 
boundary to obtain an estimate of the discontinuities B d j of 
the boundary; 

updating the set of parameters £, based upon a 
circle fitting method being applied to the continuous parts 
B ci and the discontinuities B d j of the boundary; 

repeating said step of estimating the boundary 
until the set of parameters S converges; and 

segmenting the nodule, based on estimates of the 
continuous parts of the boundary and the discontinuities of 
the boundary . 

9. The method of claim 8, wherein said constructing 
step comprises the step of increasing a radius of the 
deformable circle until the radius contacts high gradient 
points in the 3D CT volume dataset . 

10. The method of claim 9, further comprising the 
step of pre-processing a region-of -interest that 
encompasses the user- selected point using an Expectation 
Maximization (EM) based method, so as to classify and 
remove a calcification from the region-of -interest . 
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11. The method of claim 10, wherein said pre- 
processing step removes the high gradient points that 
result from the calcification of the nodule. 

12. A program storage device readable by machine, 
tangibly embodying a program of instructions executable by 
the machine to perform method steps for automatically 
segmenting lung nodules in a three-dimensional (3D) 
Computed Tomography (CT) volume dataset, said method steps 
comprising : 

receiving an input corresponding to a user- selected 
point near a boundary of a nodule; 

constructing a model of the nodule from the user- 
selected point, the model being a deformable circle having 
a set of parameters 6 that represent a shape of the nodule; 

estimating continuous parts of the boundary and 
discontinuities of the boundary until the set of parameters 
S converges, using dynamic programming and Expectation 
Maximization (EM) ; and 

segmenting the nodule, based on estimates of the 
continuous parts of the boundary and the discontinuities of 
the boundary. 
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13. The program storage device of claim 12, wherein 
the set of parameters S = [0, s] T , 0 being a position of the 
model, s being a scale of the model, and T being a 
transpose of a vector corresponding to the position 0 and 

5 the scale s of the model . 

14. The program storage device of claim 13, further 
comprising the step of representing the boundary as a sum 
B, where B= B ci ) U(U; B dj ) , B ci represents continuous parts of 

10 the boundary and B d j represents discontinuities of the 
boundary . 

15. The program storage device of claim 14, wherein 
said estimating step includes the steps of; 

15 estimating the continuous parts B c ± of the boundary 

based on a Maximum A-posteriori (MAP) estimate according to 

an equation B ci = arg maxB ci p(B C i^J, E) , I being a slice from 
the 3D CT volume dataset ; 

estimating a MAP density as p{B C i^I, E) = 1/z 
20 exp (E E (B ci ) ) , E E (B c i) being a sum of internal shape and 
external image energies, and z being a normalization 
constant ; 



minimizing the sum of internal shape and external 
image energies Eg ( B c ± ) using a time-delayed discrete dynamic 
programming method; 

connecting the continuous parts B ci of the boundary to 
5 obtain an estimate of the discontinuities B d j of the 
boundary; and 

updating the set of parameters S, based upon a circle 
fitting method being applied to the continuous parts B C i and 
the discontinuities B dj of the boundary. 

10 

16. The program storage device of claim 12, wherein 
said constructing step comprises the step of increasing a 
radius of the deformable circle until the radius contacts 
high gradient points in the 3D CT volume dataset . 

15 

17. The program storage device of claim 16, further 
comprising the step of pre-processing a region-of -interest 
that encompasses the user-selected point using an 
Expectation Maximization (EM) based method, so as to 

20 classify and remove a calcification from the region-of- 
interest . 
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18. The program storage device of claim 17, wherein 
said pre-processing step removes the high gradient points 
that result from the calcification of the nodule. 
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